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Abstract. The complex geometry underlying the Schrodinger dynamics of coherent 
states for non-Hermitian Hamiltonians is investigated. In particular two seemingly 
contradictory approaches are compared: (i) a complex WKB formalism, for which the 
r-j 1 centres of coherent states naturally evolve along complex trajectories, which leads to a 

O j class of complexified coherent states; (ii) the investigation of the dynamical equations 

for the real expectation values of position and momentum, for which an Ehrenfest 
theorem has been derived in a previous paper, yielding real but non-Hamiltonian 
classical dynamics on phase space for the real centres of coherent states. Both 
approaches become exact for quadratic Hamiltonians. The apparent contradiction is 
resolved building on an observation by Huber, Heller and Littlejohn, that complexified 
coherent states are equivalent if their centres lie on a specific complex Lagrangian 
1 manifold. A rich underlying complex symplectic geometry is unravelled. In particular 

t^- ■ a natural complex structure is identified that defines a projection from complex to real 

phase space, mapping complexified coherent states to their real equivalents. 



PACS numbers: 03.65Sq, 02.40Tt 



^ ■ 1. Introduction 

5-h ' 



We analyse the geometric structure related to complexified coherent states, that is 
Gaussian states with a formal complex centre. These states appear naturally in 
situations where the classical Hamiltonian function is complex valued, or in classically 
forbidden regions in the description of tunneling processes. Here we focus in particular 
on the quantum counterpart of complex Hamiltonians, that is the quantum time 
dependence generated by the Schrodinger equation with a non-Hermitian Hamilton 
operator. Such operators are of interest in many areas in science, in particular in physics 
and chemistry. They appear, e.g., in the description of decay processes in quantum 
mechanics, from early models in nuclear physics to the use of complex scaling in the 
computation of resonances [1] . In optics they naturally appear in the study of absorbing 
or optical active materials [2j[3], and in chemistry absorbing complex potentials are 
frequently used for numerical simulations |4j. From a more mathematical perspective 
the spectral theory of non-Hermitian operators has received renewed interest recently, 
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in particular due to questions arising from numerical analysis, such as the concept of 
the pseudo spectrum [5]. Further, the special class of non-Hermitian PT symmetric 
operators has received much attention recently, since these operators often possess a 
purely real spectrum, and have been suggested as a generalisation for the description of 
closed quantum systems (6j[T| . 

Coherent states play a crucial role in the correspondence of quantum and classical 
systems. They can be used to quantise classical systems, but also the semiclassical limit 
of a quantum system can be conveniently expressed with the help of coherent states 
[8l-fl3] . In the usual formulation quantum to classical correspondence will associate with 
a non-Hermitian Hamilton operator a Hamilton function which is complex valued, and 
therefore the corresponding Hamiltonian dynamics will generate complex trajectories in 
a complexified phase space [6| IT4] - l20] . Considering the semiclassical limit in the sense 
of the Ehrenfest theorem for expectation values, on the other hand, leads to a classical 
dynamics on real phase space developed in [21]. Here we will show that the seeming 
contradiction between these two formulations is related to an ambiguity in the definition 
of complexified coherent states. It is resolved by the identification of an equivalence class 
of coherent states with complex centres, [13j[22], beautifully encoded in the concept of 
Lagrangian manifolds. 

The plan of the paper is as follows. In the next section we summarise some of the 
complex symplectic geometry which underlies the properties of complexified coherent 
states. This will provide the framework for the following section where we consider 
quantum and classical dynamics. We focus on the discussion of quadratic non-Hermitian 
operators for which the quantum classical correspondence is exact and which form the 
basis for semiclassical considerations of more general situations. We will in particular 
compare the structures which emerge from an extension of WKB theory to the non- 
Hermitian case with an extension of the Ehrenfest Theorem [21] . 

2. Coherent states and complex structures 

It is well known that a manifold of coherent states can be interpreted as the phase space 
of the corresponding classical system [TTl[23l[2l], and how the symplectic structure of 
classical mechanics naturally arises from the geometry of coherent states. What is 
perhaps less appreciated, is that the coherent state manifold is further equipped with a 
metric and a complex structure, which is as well inherited to the classical system. As the 
metric structure does not appear in classical Hamiltonian equations of motion, it can be 
easily overlooked. This is different in the context of dissipative classical systems, where 
in addition to the symplectic flow of Hamiltonian dynamics, a metric gradient flow often 
appears. These types of dynamics are sometimes referred to as metriplectic flows [25J. 
It has recently been pointed out, how similar structures arise in the semiclassical limit 
of non-Hermitian quantum theories, where the metric of the classical phase space is 
provided by the metric on the space of coherent states [2X11261 [27] . 

Let us now recall how certain classes of Gaussian coherent states endow classical 
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phase space with a metric and a complex structure. We begin with a brief review on 
the familiar case of real coherent states. Consider a family of Gaussian states 

B = (det Im B) 1/4 ^ P . {x _ Q)+ ^ (x _ Q) . B{x _ Q)] , . 

vz\ ) (vr^)«/4 ' v > 

with Z = (P,Q) G R™ x R ra , and B e M n (C), where B is symmetric and has positive 
imaginary part, Imi? > 0. This last condition ensures that the state is in L 2 (R n ) and 
the prefactor is chosen such that the state is normalised to one. 

This coherent state manifold can in the semiclassical limit be identified with the 
classical phase space via the centre Z, and the matrix B defines a metric and a complex 
structure on phase space. The metric emerges in a natural way in a phase space 
formulation of quantum mechanics, using for example the Wigner function, see [12] 
for the following. The Wigner function of the state (CQ) is a Gaussian centred around 
Z = z, and localised on the order of H: 



H ''M = (^F <H<2 '" ^, ' G<2 '" ^, • (2 » 

where z' = (p', q') denotes the coordinate and momentum variables, and the positive 
symmetric matrix G is related to B via 



«= ' " 1 „' " r ■ (3) 





Hence the matrix G defines a metric on phase space. This metric G has the additional 
property that it is symplectic, i.e., it satisfies GQG = Q, where Q denotes the 
antisymmetric matrix 

Since QQ = —I this implies that — QGQ = G^ 1 and using this it is easy to see that 

J := —QG , (5) 

defines a complex structure on phase space, i.e., it satisfies J 2 = —I. Recall that a 
general f2-compatible complex structure on phase space is a symplectic matrix J such 
that J 2 = —I, and the matrix QJ is positive definite. For later use we note that by (J3j) 
the complex structure J can be expressed in terms of B as 



J 



ReBllmB}- 1 Im B + Re B[Im B}- 1 Re B 
-[ImB]- 1 [ImBj^ReB 




(6) 

ReB 
I 
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The physical meanings of the centre Z and the metric G become apparent when 
considering expectation values and variances of physical observables A. Let A be the 
Weyl quantisation of a smooth classical observable A, then the expectation value and 
variance of A in the state are given by 



(i/),Ai/>) _ A(r7 ^ , fA ^ 2 _ h 

2 



A(Z) + 0(h) , (AA)l = -VA{Z) ■ G' l VA{Z) + 0(h 2 ) , (7) 

i.e., Z is the centre of the phase space distribution of ip and G determines its variance. 
Thus, in the limit of h — > each coherent state collapses to a phase space point Z, and 
the matrix G encodes a local metric at this point. 

Let us now extend the previous considerations to the case that the coherent state 
is formally centred at a complex phase space point, i.e., we consider Gaussian coherent 
states on R n similar to (CQ), but with a complex centre z = (p, q) G C ra x C n , 

tf( x ) = (det Im B)V* mx _ q) ^ {x _ gyB(x _ q)] 

and B e M n (C) is again a symmetric n x n matrix with Imi? > 0. Note that while the 
centre is formally chosen complex, the wave function can still be viewed as a function of 
a real coordinate x e R n , and the condition ImB > guarantees that it is in L 2 (R n ). 
Similar states were considered previously by Huber, Heller and Littejohn, [T3 4 22], and it 
was noted that different choices of the complex centre z can lead to the same quantum 
state. In particular, it was found that two centres, z and z', define the same quantum 
state if 

z-z' eL B :={(Bq,q);qe<C n } , (9) 

where Lb is a natural complex Lagrangian space associated with the state (jHJ) which we 
will analyze in more detail below. We will show here that this result is closely related 
to the complex structure J induced by B and can be reformulated in terms of a natural 
projection from complex phase space to real phase space defined by 

Pj(z) := Rez + Jlmz , (10) 

i.e., Pj(Rez + ilm^) = Re z + Jim z, where z G C n x C n and the real and imaginary 
parts are taken component-wise. 

The main result of this section can now be formulated as follows. 

Theorem 2.1. Let ip^ be the coherent state (jHJ) and Pj be the projection f JTQj) defined 
in terms of the complex structure (Q, then 

a = ^ {z > Pj{z)) r Pj{z) , (ii) 

where with z = (p, q) and Pj(z) = (P, Q) we have 

o{z,Pj{z)) = l -{P + p).{Q-q) . (12) 
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a -21ma(z,Pj(z))/h 

W ( z >\ = ? — e -i(,'-p J (,)).G(,'-p J (, ) ) (13) 

In other words, the complex "centre" z = Re z + i Im z of the state (jSJ) is projected 
to the real centre Z = Rez + Jlmz. Hence the coherent state centred at z is physically 
equivalent to the one centred at Z = Pj(z). 

Proof. We can write a coherent state (jSJ) in the form ij){x) = Ce^ s ^ with 

S(x) = p ■ (x - q) + ~(x - q) ■ B(x - q) , (14) 

and some constant C. The crucial step is to note that this state is concentrated 
around the point where the imaginary part of S(x) is minimal, but since the parameter 
z = (p, q) e C ra x C n can be complex the minimum need not be located at x — q. Let 
us introduce Z = (P, Q) by the conditions 

V Im S(Q) = and P = V Re S(Q) = VS(Q) (15) 

then Q is the minimum of the imaginary part of S and by expanding S(x) up to second 
order around x = Q we can rewrite S(x) as 

S(x) = S(Q) + P ■ (x - Q) + ±(x - Q) ■ B(x - Q) . (16) 

The complex structure will now appear if we express Z in terms of z. We find 
VlmS'(x) = lm[p + B(x — q)] = Imp + lmB(x — Reg) — lie B Iraq and thus the 
condition VlmS(Q) = gives 

Q = Req+ [ImB^ReBlmq - [ImB^Hmp. (17) 

Since VRe5'(a;) = Re[p + B(x — q)] = Rep + ReB(x — Reg) + ImB Iraq we obtain 
further 

P = Rep-Re J B[Im J B]" 1 Imp+(Re J B[Im J B]" 1 Re J B + Im J B)Img . (18) 

These two equations yield Z = (P,Q) = Pj(z), with J given by (j5J), and hence with 
(USD we find 

= e iS ^f j{z) . (19) 

It remains to compute S(Q) — p ■ (Q — q) + |(Q — q) ■ B(Q — q) . From (fTTI) and (ITS|) we 
find B(Q — q) = P — p and hence 

S(Q) = \{P +p)-{Q-q) = <t(z, Pj(z)) . (20) 
The form of the Wigner function (fl3l) follows from (flT!) and (J2]). □ 
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In the remainder of this section we want to elucidate the complex symplectic 
geometry underlying and connecting the complex matrix B in the definition of a 
coherent state ( ITT]) , the complex structure J (EJ), and the Lagrangian submanifold Lb 
(J9]). Obviously J and Lb are both defined in terms of B. We can further show that 
there are one-to-one relationships between all three of them. 

Let us recall that a linear subspace L C C n x <C n is called Lagrangian if Q\l = 
and dim L = n, and positive Lagrangian if in addition the quadratic form 

h( z , z') := ]-z ■ Qz' (21) 

is positive on L, i.e., h(z, z) > for all z G L. It is a well known result [28] that any 
positive Lagrangian subspace can be written in the form fl9]): 

Lemma 2.2. The subspace Lb = {(Bq,q); q G C n } defined in is Lagrangian if B is 
symmetric, and positive Lagrangian iflmB > 0. On the other hand, i/LcC"xC™ is 
a positive Lagrangian subspace then there exists a symmetric B G M n (C) with ImB > 
such that L = {(Bx,x) ; x G C d } . 

Proof. It is clear from the definition that dimL = n. To check that Q\l = we choose 
z = (Bx, x) G L and z' = (Bx', x') G L and find z ■ Qz' = —Bx ■ x' + x ■ Bx' = x ■ (B T — 
B)x' = 0, since B is symmetric. To check positivity we consider hi(z,z) = iz ■ Qz/2 
with z = (Bx, x) which gives 

h L (z,z) = l -[-{Bx) -x + x ■ Bx} = l -x ■ [B - B T ]x = x-lmBx > . (22) 

Now assume L to be a positive Lagrangian subspace and consider the projection 
7r : L — )• C n defined by 7i(p,q) = q. Then ker7r = {0} because if z = (p,q) G ker7r, 
then q = and hence iz ■ Qz/2 = 0, thus positivity of L implies z = 0. Therefore the 
map 7r is invertible and since it leaves the q component invariant the inverse must be of 
the form 7c~ 1 (q) = (Bq,q) for some matrix B, i.e., L = {(Bq,q) ; q G C n }. That B is 
symmetric and has positive imaginary part follows now as before from the fact that L 
is positive and Lagrangian. □ 

This establishes the one-to-one correspondence between Lagrangian subspaces and 
complex symmetric matrices with positive imaginary part. Let us now relate complex 
structures and positive Lagrangian subspaces. By ( TTOl) and (flTil . the complex centres 
z and z' define the same state if Pj(z — z') = 0. Hence, the set of equivalent complex 
centres is given by 

L := ker Pj = {z G C n x C n : Pj{z) = } . (23) 

According to the work of Heller et. al. JT3l[22] we expect that L = Lb- Let us, however, 
first show that L is actually a positive Lagrangian manifold, and furthermore, that the 
set of ^-compatible complex structures is isomorphic to the set of positive Lagrangian 
subspaces of C x C". 
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Lemma 2.3. Let J be a Q- compatible complex structure on R n x R™ (see the definition 
after equation (jSJ) ) and define 

Pj(z) :=Rez + Jlmz , (24) 

then 

L : = ker Pj = {z G V c ; Rez + Jlmz = 0} (25) 

a positive Lagrangian subspace. Conversely, for every positive Lagrangian subspace L 
there exists a compatible complex structure Jl such that L = ker Pj L , i.e., 

z G L <£> Rez + J L lmz = . (26) 

Proof. Note that since J 2 = — I the relation Re 2 + Jlmz = can be rewritten as 

lmz = JRez, (27) 

i.e., z G L means 2 = (J + iJ)Rez. Since J is non-degenerate we clearly have 
dim c L = n, and for z,z' G L we get 

* • fb' = Re z ■ (I + U T )tt(I + i J) Re z 

= Re z ■ (SI - J T QJ) Re z' + i Re z ■ (J T Q + ft J) Re / 

and if J is symplectic and f2J = G symmetric we get that z ■ flz' = 0, and thus L is 
Lagrangian. Furthermore we find for z — (I + iJ) Re z G L 

^z ■ Qz = Rez ■ GRez , (29) 

hence L is positive. 

On the other hand, assume L C C n x C n to be a positive Lagrangian subspace 
and consider the map Re^ : L — > R n x R n , defined by Rei(z) = Rez. We claim 
that this map is invertible. To see this assume z G kerRe^, i.e, Re z = 0, then 
\z ■ Qz = ilm z ■ film z = 0, hence z = by the positivity of L, so kerRe^ = {0} 
and Re^ is invertible as claimed. The inverse must be of the form Re^ 1 (f) = v + iJv 
for a linear map J : R" x 1" -> R" x R" Then (j2HD with Re 2; = v shows that if L 
is Lagrangian J must be symplectic and G := flJ symmetric, and (|29|) shows that G 
must be positive. Then J 2 = QGQG = Qfl = —I, therefore J is a compatible complex 
structure. □ 

In summary, we have shown that the set of complex symmetric matrices with 
positive imaginary part, the set of positive Lagrangian subspaces, and the set of Q- 
compatible complex structures are all isomorphic to each other. What we have not 
shown yet is that Lb is actually mapped to the complex structure (jEJ), i.e., that 



ker Pj = L B 



(30) 
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Since dim ker Pj = dim Lb it is enough to show that Lb C ker Pj, i.e., that for any 
z G Lb we have Rez + Jlmz = 0. Now any element in Lb is of the form z = (Bq, q) 
for some q G C ra and a short calculation gives 





ImB ReP 
/ 




(31) 



and hence z G ker Pj for all z E Lb means 




imP\ T ImB ReB 

o +J / 



. 



(32) 



Solving this equation for J then gives the expression (J6]) which we have already 
encountered. Hence (I30I) holds. 

For completeness we finally note that the metric G = QJ defines a Kahler structure 
on complex phase space which turns Pj into an orthogonal projection: 

Lemma 2.4. Let h(z, z') := z ■ Gz' — \z ■ Qz' be the hermitian inner product on C n x C n 
defined by G = QJ, then Pj is the unique projection onto R n x R n which is hermitian 
with respect to h, i.e., h(Pjz,z') = h(z,Pjz'). 

Proof. Pj is a projection, so it is hermitian with respect to h(z, z') if the kernel and 
image are orthogonal to each other. Then it is as well uniquely determined by its 
image. Since by ( 1271) any z G L = ker Pj is of the form z = (I + iJ)x for some 
x G R n x R™, we get for z = (I + iJ) G L and z' = x' G R" x R n = ImPj that 
h(z,z') = x - (I + U)\G - iQ)x' = x ■ [G + J l Vt + i(J*G - Q)]x'. But since G is 
symmetric G = QJ implies G = —J l Q and from GJ = —Q we obtain J l G = Q, 
therefore h(z, z') = for all z G ker Pj and z' G Im Pj. □ 

We have shown that coherent states with a complex centre are organised along 
Lagrangian submanifolds of physically equivalent coherent states one of which has a 
real centre. In what follows we shall investigate the time dependence of these structures 
under the evolution with non-Hermitian Hamiltonians. In particular, we will focus 
on the analytically solvable case of quadratic Hamiltonians, which lies at the heart of 
semiclassical considerations for more general systems. 



3. Schrodinger dynamics with complex quadratic Hamiltonians 

Here we will investigate the Schrodinger dynamics generated by complex quadratic 
Hamiltonians that are given as Weyl quantisations of complex quadratic forms on 
phase space. For these Hamiltonians semiclassical approximations are exact, and we 
restrict ourselves to these purely quadratic Hamiltonians to understand the essence of 
the dynamics in detail. It is straightforward to include also linear terms; here, however, 
we want to keep the discussion concise. 
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We continue to denote by z = (p, ?) e R" x R ra points in phase space and set 

U{z) = ^z- Hz , (33) 

where H 6 Af 2n (C) is a complex symmetric 2n x 2n matrix and the quantum 
Hamiltonians we will consider are given by the Weyl quantisation of quadratic functions 
of the form H, 

h 1 ih 

^pp^x + -rx ■ H qp V x + -x ■ HggX - — tr H qp , (34) 







H = 


- h -v 

2 


( Hpp 


H pq \ 


\H qp 


Hgg) 



where H = J . We will in general allow the matrix H to be time dependent 

yflgp tlgg J 

without explicitly indicating this in the notation. Our aim is to study the solutions to 
the time dependent Schrodinger equation 

ihd t ip = Hip , (35) 

for initial states given by coherent states. Since our Hamilton operator is in general 
not self-adjoint the question of whether this equation has solutions in suitable function 
spaces is not trivial. To illustrate the issue, consider the following simple example: If 
the Hamiltonian is given by H{z) = iq 2 /2 the time evolution operator is of the form 
U (t) = e"^ x2 and taking for instance an initial state of the form ipo(x) = e~^ x2 it follows 
that 

ip(t } x) =e~^ x2 (36) 

and hence ip(t,x) £ L 2 (R) for t > b. 

Problems of this kind are avoided if the imaginary part of H is chosen to be non- 
positive. For ImH < the Schrodinger equation generates a contracting semigroup, and 
the quadratic case has been studied in some detail. We mention [29], where the Weyl 
symbols of the time evolution operator have been constructed explicitly using complex 
symplectic geometry, and |30j for some early rigorous results on the damped harmonic 
oscillator. Here we will further analyse the consequences on the geometric structures we 
have highlighted in the previous section. In addition, the case of non-negative ImH is 
often of interest, in particular in the context of PT-symmetric quantum systems. Thus, 
we allow for general complex H here, but we only consider special initial conditions for 
which explicit solutions can be computed, at least for short times. 

We will investigate the dynamical behaviour of initially Gaussian coherent states 
that is generated by a Hamiltonian operator of the form (|34"|) . Similar to the real valued 
case, the class of Gaussian coherent states, now with a complex centre, is invariant under 
this time evolution, as we shall see in the following. For this purpose we consider time 
dependent Gaussian coherent states of the form 

^ t ,x) = e te (*)^^^^eHW*)-(*-^))+i( a! -^))- fl ( t )( !B -9W)] = e^Vf^Or) , (37) 
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where z(t) = (p(t),q(t)) G C n x C n , B(t) G M n (C) is symmetric and has positive 
imaginary part, Im B(t) > 0, and a(t) G C Inserting the state fl37|) as an ansatz into 
the Schrodinger equation ( 1351) and separating terms with different powers of (x — q) 
yields the following set of differential equations for (p(t) , q(t)) , B(t) and a(t): 

-p + Bq = H' q + B% (38) 

-b = % q + n; q B + b% v + Bn; p B (39) 

-d + ~ tr OBiT 1 ) = -l\p. q -U\- -[tvU^ + tr(?£B)] , (40) 

where denote derivatives of 7i(z) with respect to p, and p and q, etc.. If 

we choose p and g to be solutions to Hamiltons equations, i.e., p = —7~L q and q = H' , 
then the first equation is satisfied, and furthermore using the second equation we can 
simplify the third, thus arriving at the simplified system 

z = VLHz (41) 
B = —H qq — H pq B — BH qp — BH PP B (42) 

a = I[p - q - U{z)} + i tr[H pp B - H qq B^ 1 } . (43) 

Here the first equation is Hamilton's equation with a complex Hamilton function and 
the third equation can be integrated once the first and the second are solved. The 
solutions to the second equation can be obtained most easily using symplectic geometry 
which will be reviewed in what follows. This set of equations is a complex extension of 
the classical approach to coherent state propagation of Hepp, [9j, and Heller [ID] , which 
is used and developed further in many areas (see, e.g., the review [12] or [31] for an 
overview of more recent mathematical developments). 

For complex H equation (141j) leads to complex solutions z(t), even if the initial 
condition is chosen to be real, and thus we will obtain coherent states with complex 
centres. As discussed in the previous section a complex centre has no direct physical 
meaning, but using a complex structure it can be projected to a physically meaningful 
real centre. We will now apply the complex symplectic geometry we developed in the 
last section to understand the relation between the dynamics of the complex centre and 
its projection to real space. 

In a previous paper [21J we concentrated on the dynamics of the Wigner function 
which directly yields the expectation values and hence the real centre of a state. This 
considerations led to a non-Hermitian version of Ehrenfest's theorem with a new type of 
classical dynamics emerging in the semiclassical limit. We derived an evolution equation 
for the Wigner function, which in the case of a quadratic Hamiltonian reduces to 

hd t W(t, z) = -(- ^-A ImH - hz - Re HttV -2z-lm Hz^j W(t, z) , (44) 

where all derivatives are with respect to z, and 

A lmH := -V-fi T Im£mV . (45) 
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While the evolution equation in [2T] for general Hamiltonians is a semiclassical 
approximation, the quadratic case flH|) is exact. 

If ip is of the type ( 1371) the Wigner function is of the form 

W u z ) = 2_^ e -fa-z<t))-G(t)(*-m) (46) 
{ r nn) n 

with Z(i)el"xE",a symmetric G(t) G M 2n (R), and 0(t) G R. Inserting the ansatz 
(146p into equation ( I44p . and separating different powers of (z — Z), leads to the following 
set of equations 

Z = VlReHZ + G- 1 lmHZ (47) 
G = ReH{lG -GttReH -ImH + Gtt T ImHttG (48) 

$ = --Z-ImHZ- -tr[lmHQGQ T ] (49) 
1% — 

It can be verified, that this set of equations is also compatible with the dynamical 
equations f )38|) . f )39|) . and PU|) obtained from the coherent state ansatz in the Schrodinger 
equation, if we demand p and q to be real. Thus, equations ( l4il . ( )42|) . and ( l4Uj) are not 
the unique dynamical equations for the propagation of coherent states for non-Hermitian 
Hamiltonians. 

The two different sets of equations that we have obtained, f|4Tl) . (j4~2l . and f|47j) . 
fj48|) . are supposed to describe the dynamics of the same physical state. In what follows 
we will discuss how they can be related using complex structure associated with the 
coherent states. 



3.1. Symplectic evolution 

To solve the evolution equations obtained above, in particular the nonlinear matrix 
Ricatti equations (142]) and ( l4~8l) . we have to understand how the geometric structures 
discussed in the previous section evolve in time under the action of complex Hamiltonian 
dynamics. For this purpose, we first investigate the action of a linear symplectic map 
on a positive Lagrangian subspace L, i.e., we change L to SL with S G Sp(n, C). 
Here Sp(n, R) and Sp(n, C) denote the set of real or complex 2n x 2n matrices S with 
S T QS = Q, i.e., the real and complex linear symplectic groups. Since any z G SL is of 
the form z = Sz Q for some Zo G L we get z ■ VLz' = z ■ S T QSz' = zo ■ £lz' = 0, since L is 
Lagrangian, and thus SL is, too. Furthermore 

l -z' ■ Qz' = l -z ■ S T QSz , (50) 

thus, if 5* = S, i.e., S G Sp(n, R), then SL is positive, too. If S is complex, SL does 
not have to be positive any more. 

We will mainly consider situations in which S is the solution to Hamilton's equation, 
i.e, S(t) satisfies 

S = QHS , with S(t = 0) = / , (51) 
where H G M2 n (C) is symmetric. 
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Lemma 3.1. Assume L to be a positive Lagrangian subspace and S(t) a solution of 
(J5l]) . Then there exists a T HL such that for all t G [0,T HL ) S(t)L is again a positive 
Lagrangian subspace. IflmH < we can take Th,l = oo. 

Proof. Since S(t) is close to the identity for small t , S(t)L will be positive by continuity 
for sufficiently small t. If ImH < we proceed as follows. With Sz G SL for z G L we 
have to consider i(S'z) • QSz/2 = iz ■ S T QSz/2 for z G L. From ( )5T|) we find 

^U^OSJ = ^S T [-HQn + QQH]S = -S T lmHS. (52) 

Thus, if ImH < then ^z-S^fiSz > and S(t)L is therefore positive for all t > 0. □ 

Let us now investigate how S and the complex structure transform if we apply a 
symplectic map to L. 

Proposition 3.2. Let L be a positive Lagrangian subspace and S G Sp(n, C) such that 
SL is still positive. Then 



B SL = S*B L (53) 



where the action of S = ( *^ P9 ) on B L is defined by 

\Sqp S qq ^ 

S*B L := (S pp B L + S pq )(S qp B L + Sgg)^ 1 , (54) 

(ii) and 

J SL = (ReS -lmSJ L )J L (ReS -IiiiSJl)- 1 (55) 
G SL = ^(Re S - Im SJ L )ft T G L (Re 5 - Im SJl)" 1 . (56) 

Proof. Let z G L, then there exists a q G C n such that z = (BLq,q), by Lemma 
12. 2[ and since Sz G SX there exists a q' G C™ such that Sz = (B SL q',q'). Now 
Sz = S(B L q, q) = (S pp B L q + S pq q, S qp B L q + S qq q) and hence we obtain the two equations 

{S P pB L + S pq )q = B S Lq' , (S qp B L + S qq )q = q' . (57) 

From the second equation we get q = (S qp B L + S g? ) _1 g' and inserting this into the first 
gives (S pp Bl + Spq)(S qp BL + S qq )~ l q' = BsLq', which is the first result. 

To derive the second result we note that z G L means z = Re z + iJ^ Re z, by f j27|) . 
and similarly Sz G SL means Sz = Re(Sz) + iJsLRe(Sz) and thus we arrive at the 
expressions 

Im[Sz] = J SL Re[Sz] = J SL (Re S - Im S J L ) Re z (58) 
Im[Sz] = Im[S(Re z + i J L Re z)] = (Im S + Re S J L ) Re z . (59) 

Comparing these two expressions for Im(Sz) gives Jsl = (ImS + Re SJl)(Rg S — 
ImSJ^) -1 and with Jf = —1 we furthermore obtain (ImS + ReSJ^) = (ReS — 
Im SJl)Jl- The result for Gsl then follows from Gsl = —QJsl- □ 
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Note that there is a certain similarity in the structure of equations (I53p . (|55p and 
f l56|) . In fact we can rewrite (155]) and fl56l) as 

J SL = (ReSJ L + lmS)(-lmSJ L + ReS)- 1 (60) 

G 5 l = {ftReSQ T G L + QlmS){-lmSQ T G L + ReS)- 1 (61) 



= ®*Jl and G 5L = $*G L , (62) 



, ReS lmS\ , ^ /^ReSlF filmS . 

*=' ImS Re J ^ $= -Im^ ReS ' " (63) 



then 
with 



If the symplectic matrix S is a solution of the differential equation (1511 then this 
induces corresponding differential equations for the evolution of the matrices B$l and 
Jsl which we shall now derive. 

(h H \ 

Theorem 3.3. Let S(t) be a solution to ( l5Tj) with H = pp P9 ) , and L a positive 

yflqp tl qq J 

Lagrangian subspace, then there exists a Th,l > such that S(t)L is positive for 
t G [0,T HtL ] and we have 

Bsl = —HqpB S L — B SL H pq — H qq — B SL H pp B S L (64) 

(ii) and 

j SL = n Re H J SL - J SL tt Re H + tt Im H + J SL ^ Im HJ SL (65) 

G SL = Re HVlGsi - G SL ^ ReH-lmH + G SL ^ T Im HttG SL (66) 

Furthermore iflmH < we can take Th,l = oo. 

Proof. Since S(0) — I it is clear that for small t the space S(t)L will still be positive, 
hence there exists a T H L such that SL is positive for t e [0, T HjL ] . Now from (l5Tj) we 
get 

Spp ^pq \ | H q pSpp H qq Sqp H qp S pq H qq S qq 

Sqp ^qq J V HppS pp ~\~ H pq S q p HppSpq ~\~ HpqSqq 

then differentiating the relation (|53|) and using (|67|) gives 

-Bsl = (S pp Bl + S pq )(S qp BL + 5 gg ) 1 — BsL(S qp BL + S qq )(S qp BL + S^) 1 

= ~Hq p (SppBx J + Spqj^SgpBxj + S'gq) — Hqq^SqpBj^ + Sqqj^SqpB]^ + iSgg) 

— B SL H p p(SppB L + S p q)(S qp B L + S qq ) 1 — B SL H pq (S gp B L + Sqq)(S qp B L + Sgg) 1 

= —H q pBsL — H qq — BsiHppBsL — B S hB vq 

(68) 



(67) 
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To prove the second set of relations we first rewrite (I55p as J$l — AJ^A^ 1 with 
A = Re S — Im SJl, then 

j SL = AJlA" 1 - AJlA^AA^ 1 = AA^Jsl - JslAA' 1 . (69) 



Then from (EI]) we find ReS = ttReH ReS - ttlmH ImS and ImS = ttlmHReS + 
QReH Im S and using these relations we find 

A = ReS -lmSJ L = [nReH -nimHJ SL }A (70) 

and this leads to 

j SL = ttRe HJ SL + ttlmH- J SL tt ReH + J SL tt Im HJ SL ■ (71 ) 

The result for Gsl then follows using the relation Gsl = ttJsi □ 



The formal similarity of the equations ( 1661) and ( |64j) suggests to define a Hamiltonian 
K((, z) on the doubled phase space by 

r _ , 1_ ,f-tt T lmHtt ttReH\f(\ 

then the matrix $(t) from (J63l) satisfies 

• = /0 -A /-^Imtfn fiReif\ 

I J y I -Rei/fi ImH I ' K J 

And so by solving (173"]) with $(£ = 0) = / we find a matrix such that 

G(t) = $(t),G (74) 

is a solution to (08} with G(t = 0) = G. 

3.2. Quantum evolution 

The results from Theorem 13.31 allow us to solve the non-linear Riccati equations (T4"2l 
and (j4*8l) in terms of solutions to linear Hamiltonian equations, which we will exploit in 
what follwos. 

We first consider the Schrodinger equation for a coherent state in position 
representation, (I3TI) . Let S{t) G Spin, C) be the solutions to 

S = ttHS , with S(0) = I (75) 

then z{t) = S(t)zo is a solution to (jUJ) and by Theorem 13. 3 [ part (i), and Proposition 
13. 2\ part (i), S*B is a solution to f|42|) . Hence we conclude 
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Theorem 3.4. Let L = Lb be a positive Lagrangian subspace, then there exists a 
Th,l > such that for t e [0,Th,l) the solution to the Schrddinger equation with 
ip{t = 0) = is given by 

ik(t) = e^Vfgf (76) 

where a G C is the solution to (|4"3"|) with a(0) = 0. 

The phase factor is related to the action along S(t)z and also contains Maslov-phase 
type contributions. 

The matrix Sfy^B defines a time dependent complex structure J(t) via ([6]) which 
projects the complex centre S(t)z to the real centre 

Z{t) = P J{t) {S(t)z) (77) 

and we can use Theorem 12. II to express the Wigner function in terms of projections from 
the complex dynamics S(t). 

Alternatively we can solve the purely real set of equations (HHjl and (1471) to directly 
obtain the motion of the real centre. Let $(t) be the solution to (!73|) and Z(t) a solution 
to dHJ) with G(t) = $(t)*G then we have 

Theorem 3.5. Let G be a symplectic positive definite symmetric matrix. Then there 
exists a Th,g > such that for t e [0, Th,g) the unique solution to the Wigner von 
Neuman equation (Hlfl with initial condition W(z) = , ^ n e~H( z ~ z )' G ( z ~ z ) is given by 

W a, z ) = ^— Q -U*- z V)>m)*G]{ Z -z{t)) (78) 

[7rh) n 

where (3(t) Glisa solution to ( 149]) with (3(0) = 0. 

One of the characteristic features of the dynamical equation (jUJ for the real centre 
Z(t) is that it is in general not autonomous, the coefficients of this equation will depend 
on t via the metric G(t). However, in many cases there are special solutions for which 
the metric is time independent, corresponding to fixed points of the evolution equation 
( |48l) . To analyse the possible time independent complex structures we have to set the 
expression for the time derivative of the metric G(t) in ( T4"8]) . or equivalently the time 
derivative of the matrix B(t) in ( )42|) . to zero. Thus we obtain quadratic matrix equations 
for the fixed points Go and B , respectively. Let us illustrate this observation with a 
few examples. 

(1) Assume the Hamiltonian is anti-Hermitian, i.e., Re H = 0, then f|4*8|) with G = 
becomes ImH = GQ T ImHQG, and if we assume furthermore that for some 7 > 
we have ImH = —^S, where 5* is symplectic, symmetric and positive, then we find 
that G = S (here we used that Q T SQ = S^ 1 ). The assumptions on ImH hold for 
instance if n = 1 and ImH is negative definite (with 7 = detlmif). Thus in this 
case the metric and the associated complex structure are constant and the equation 
of motion for the centre simplifies to 

Z = -S^-fSZ = - 7 Z . (79) 
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Hence we find a uniform contraction towards the origin. For this Hamiltonian we 
can as well solve f !73|) explicitly and using f!74|) we find that for an arbitrary initial 
Go the solution to ( j48i) is given by 

G(t) = (G + tanh(7t)5)(tanh( 7 t)G + S^S = S + 0(e" 7 *) , (80) 

hence the stationary solution G = S we found above is a global attractor, to which 
any other solution converges exponentially fast to. Note as well that G(t) can be 
extended to some negative t but will eventually become singular. 

(2) We previously discussed the example T-L(z) = iq 2 /2 where the time 
evolution will be only defined for finite time. For the discussion of this case it 
is most convenient to use (H21 . which gives B = — i and hence B(t) = Bq — itl. 
Since Im B(t) = lmB — tl we see that the condition 1m B(t) > holds only for 
a finite time, after which the corresponding metric will blow up. The equation for 
the centre can also easily be solved; assume for simplicity that KeB = 0, then 
P = and the position reaches infinity in finite time 

Q(t) = j^Qo • (81) 

(3) We now have a look at a harmonic oscillator with damping induced by a momentum 
dependent imaginary part. We choose 

^(P, <?) = yP 2 + y <7 2 , (82) 

where the parameter 5 G C is assumed to satisfy \5\ = 1 and Re 5, Im<5 > 0, hence 
Im"H(p, q) = — Tie 51m 5 p 2 < 0. Therefore 5 parametrizes the strength of the 
damping relative to the kinetic energy. Note that choosing \5\ ^ 1 just amounts to 
rescaling of u — > \8\cu and t — > \5\t. Using ([421 we find that B = iu5 is a constant 
solution with Imi? > 0. We can then determine the corresponding metric G and 
the equations of motions for the centre which read 

p = — tu 2 q — 2lo Im5p , q = P- (83) 

For comparison with the classical damped harmonic oscillator we transform this set 
of first order equations into a second order equation for q, 

q + 2ulm5 q + u 2 q = . (84) 

We see that due to the metric this describes an underdamped oscillator, since 
lm.5 < \5\ = 1, irrespective of the choice for S. 

(4) It is instructive to include an example with a linear term in z = (p, g) £ R 2 , 

H 7 {z) = -z-z + i^-Qz (85) 
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where 7 G R 2 . The inclusion of Q in the linear term is convenient, it implies that 
if z is to the right of 7 the term is negative and we have damping, and if z is to the 
left of 7 the term is positive and we have enhancement. This is a PT symmetric 
system. It can be brought to the more familiar form q) = \p 2 + V(q) with 
V(—q) = V(q) via a canonical rotation of the phase space variables. Since Re H = I 
and ImH = we find that G = I is a solution for all times and with this initial 
choice the equation of motion for the centre Z, see [21], becomes Z = Q(Z — 7). 
Hence Z(t) = 7 + O(t)(Z - 7) with 0(t) G SO (2) denoting a rotation by t. We 
can as well solve the equation for (3 and find (3(t) = —7 • (Z(t) — Z ). Thus, the 
centre of the Wigner function evolves along circles as for the harmonic oscillator, 
but the circles are shifted due to damping and enhancement in different parts of 
phase space. The relation to a real harmonic oscillator "Ho can be directly seen in 
the following way. Introducting the complex translation T{^) := e~ 7 ' 2 we have 
T( r y)~ 1 'H T('y) = "H 7 — |7| 2 /2, and thus the operator is conjugated to the harmonic 
oscillator by a non-unitary operator, and thus the spectrum is purely real. The 
norm of the state stays bounded over time although it oscillates, which reflects the 
fact that the eigenvalues of the Hamiltonian are real, but the eigenf unctions are 
not orthogonal for 7 7^ 0. 

4. Summary 

Coherent states are a useful tool for the investigation of semiclassical limits of quantum 
theories. The investigations presented here can be viewed as part of a programme to 
understand the classical dynamics emerging from the semiclassical limit of general non- 
Hermitian operators. We recently formulated an Ehrenfest Theorem for non-Hermitian 
operators [2T] . in which the classical dynamics is given by a combination of a symplectic 
and a metric gradient field, which are generated by the real and imaginary part of the 
Hamilton function, respectively. This is a very different type of dynamics compared to 
what one would expect from extending standard WKB theory to complex Hamiltonians, 
which results in a Hamiltonian flow on complexified phase space. The main result here, 
is the proof that these two approaches are physically equivalent and are related by a 
projection from complexified phase space to real phase space, 

i^J, (86) 

where J is a complex structure on phase space which is determined by the physical 
states and becomes a dynamical variable in our theory. 

We restricted ourselves to quadratic Hamiltonians and Gausssian coherent states 
here, because both semiclassical approaches become exact in this case, and we could 
focus on the complex symplectic geometry relating them. This will form the basis for 
extension to more general systems following [21J. It is well known that semiclassical 
methods based on dynamics in complexified phase space often run into difficulties 
related to analytic extensions, e.g., complex trajectories often develop singularities, and 
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Hamiltonians which are very close on real phase space can have very different analytic 
extensions. The results described in the present paper and in [21] provide an alternative 
approach which is non-Hamiltonian but manifestly real and thus avoids these problems. 
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